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We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent 
magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused 
by the three kinds of the inhomogeneities, i.e., inhomogeneous turbulence; the nonuniform fluid 
density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the 
well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of 
the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in 
the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy 
and magnetic pumping effects which are due to the effect of the mean magnetic field on the large- 
scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be 
stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition 
field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context 
of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent 
magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. 
In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale 
dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic 
field, while they do not affect the turbulent magnetic diffusion of the poloidal field. 
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1. Introduction 



Magnetic fields observed in astrophysical plasma are 
strongly inhomogeneous (see, e.g., Moffatt 1978, Parker 
1979, Krause and Radler 1980, Zeldovich et al. 1983, 
Ruzmaikin et al. 1988, Stix 1989, Roberts and Soward 
1992, Kulsrud 1999, and references therein). For in- 
stance, the sunspots and the solar active regions are re- 
lated to the strongly inhomogeneous large-scale magnetic 
fields. The scales of the magnetic inhomogeneities, e.g., 
in the Sun are much smaller than the radius of the Sun 
and usually much larger than the size of granules of the 
solar convection. One of the mechanisms of the forma- 
tion of the magnetic inhomogeneities is associated with 
the magnetic buoyancy instability of stratified continu- 
ous magnetic field (see, e.g., Parker 1966, Gilman 1970, 
Priest 1982). The magnetic buoyancy instability is ex- 
cited if the scale of variations of the initial magnetic field 
is less than the density stratification length. This mecha- 
nism does not include explicitly magnetic flux tubes. On 
the other hand, the buoyancy of the magnetic flux tubes 
as a mechanism of the formation of the magnetic struc- 
tures was studied in a number of publications (see, e.g., 
Parker 1955, Spruit 1981, Spruit and van Ballegooijen 
1982, Schiissler et al. 1994, Moreno-Insertis et al. 1996). 
Note also that the problem of the storage of magnetic 
fields and the formation of flux tubes in the overshoot 
layer near the bottom of the solar convective zone was 
investigated, e.g., by Spiegel and Weiss (1980), Tobias 
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et al. (2001), Tobias and Hughes (2004), Brandenburg 
(2005). 

Another universal mechanism of the formation of the 
nonuniform distribution of magnetic flux in flows of the 
conducting fluid is a magnetic flux expulsion. In particu- 
lar, the expulsion of magnetic flux from two-dimensional 
flows (a single vortex and a grid of vortices) was demon- 
strated by Weiss (1966). In the context of solar and stel- 
lar convection, the topological asymmetry of stationary 
thermal convection plays very important role in the mag- 
netic field dynamics (Drobyshevski and Yuferev 1974). 
Fluid rises at the centers of the convective cells and falls 
at their peripheries. This results in that the ascend- 
ing fluid elements (contrary to the descending fluid ele- 
ments) are disconnected from one another. This causes a 
topological magnetic pumping effect allowing downward 
transport of the mean horizontal magnetic field to the 
bottom of a cell but impeding its upward return (Droby- 
shevski and Yuferev 1974, Zeldovich et al. 1983, Gal- 
loway and Proctor 1983). The fine structure of a sunspot 
is determined by the local interaction between magnetic 
fields and turbulent convection near the Sun's surface. It 
was shown recently by Thomas et al. (2002) that a down- 
ward pumping of magnetic flux may cause of filamentary 
structures in sunspot penumbrae. In particular, the mag- 
netic field lines are kept submerged outside the spot by 
turbulent, compressible convection, which is dominated 
by strong, coherent, descending plumes. 

Turbulence causes additional effects, e.g., the turbu- 
lent diamagnetic and paramagnetic drift velocities of 
the mean magnetic field (Zeldovich 1956, Krause and 
Radler 1980, Vainshtein and Kichatinov 1983, Kichati- 
nov 1991, Kichatinov and Rudiger 1992, Kichatinov and 
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Pipin 1993, Kleeorin and Rogachevskii 2003, Radler et 
al. 2003, Rogachevskii and Kleeorin 2004). In particu- 
lar, an inhomogeneity of the velocity fluctuations leads 
to a transport of mean magnetic flux from regions with 
high intensity of the velocity fluctuations (turbulent dia- 
magnetism, see, e.g., Zeldovich 1956, Krause and Radler 
1980, Vainshtein and Kichatinov 1983, Kichatinov and 
Riidiger 1992, Radler et al. 2003). On the other hand, 
an inhomogeneity of magnetic fluctuations due to the 
small-scale dynamo causes turbulent paramagnetic veloc- 
ity, i.e., the magnetic flux is pushed into regions with high 
intensity of the magnetic fluctuations (Vainshtein and 
Kichatinov 1983, Kichatinov 1991, Radler et al. 2003). 
Another effects are the effective drift velocities of the 
mean magnetic field caused by inhomogeneities of the 
fluid density (Kichatinov 1991, Kichatinov and Riidiger 
1992) and pressure (Kichatinov and Pipin 1993). In a 
nonlinear stage of the magnetic field evolution, inhomo- 
geneities of the mean magnetic field contribute to the 
diamagnetic or paramagnetic drift velocities depending 
on the level of magnetic fluctuations due to the small- 
scale dynamo and level of the mean magnetic field (Ro- 
gachevskii and Kleeorin 2004) . The diamagnetic velocity 
causes a drift of the magnetic field components from the 
regions with a high intensity of the mean magnetic field. 

The pumping of magnetic flux in three-dimensional 
compressible magnetoconvection has been studied in di- 
rect numerical simulations by Ossendrijver et al. (2002) 
(see also review by Ossendrijver 2003). The resulting 
magnetic pumping effects are isolated in the direct nu- 
merical simulations by calculating the turbulent diamag- 
netic and paramagnetic velocities. The pumping effect in 
the vertical direction is found as a predominating down- 
ward advection with a maximum speed in the turbulent 
convection of about 10 percents of the turbulent velocity 
(Ossendrijver et al. 2002). 

The turbulent diamagnetic and paramagnetic veloci- 
ties were determined analytically in previous studies only 
for purely hydrodynamic turbulence. A relation to the 
turbulent convection was made in some studies (see, e.g., 
Kichatinov 1991, Kichatinov and Pipin 1993) only phe- 
nomenologically, using the equation (u' 2 ) cx gTo{u' z s') 
which follows from the mixing-length theory. Here (u' z s') 
is the vertical turbulent heat flux, u' and s' are fluctua- 
tions of fluid velocity and entropy, g is the acceleration 
of gravity and To is the characteristic correlation time of 
turbulent velocity field. This relationship implies that 
the vertical turbulent heat flux plays a role of a stirring 
force for the turbulence. However, a more sophisticated 
approach implies a solution of a coupled system of dy- 
namical equations which includes the equations for the 
Reynolds stresses (u^u'j), the turbulent heat flux (s'u-), 
the entropy fluctuations (s' s'), the magnetic fluctuations 
(bi bj), the cross helicity tensor (bi u'j) and (bi s') in a tur- 
bulent convection. The latter has not been taken into ac- 
count in the previous studies of the small-scale magnetic 
buoyancy and magnetic pumping effects caused by the 
turbulent diamagnetic and paramagnetic drift velocities. 



Note that the turbulent convection can strongly affect 
these phenomena. 

The goal of this study is to determine the nonlin- 
ear drift velocities of the mean magnetic field in a tur- 
bulent convection. We demonstrate that the nonlinear 
drift velocities depend on the different kinds of the in- 
homogeneities: (i) inhomogeneous turbulence; (ii) the 
nonuniform fluid density and (iii) the nonuniform tur- 
bulent heat flux. The inhomogeneous turbulence causes 
the well-known turbulent diamagnetic and paramagnetic 
velocities. In addition, the nonlinear drift velocities re- 
sults in the small-scale magnetic buoyancy and magnetic 
pumping in the turbulent convection. These phenom- 
ena are different from the large-scale magnetic buoyancy 
and magnetic pumping effects. The large-scale phenom- 
ena are caused by the influence of the mean magnetic 
field on the large-scale fluid flow. Our study shows that 
these large-scale effects are stronger than the small-scale 
magnetic buoyancy and magnetic pumping only for a 
strong mean magnetic field (about equipartition field). 
We study the small-scale magnetic buoyancy and mag- 
netic pumping effects in the context of the solar and 
stellar turbulent convection. In particular, we demon- 
strate that in the main part of the solar convective zone 
the small-scale magnetic pumping effect dominates, while 
near the solar surface the radial drift velocity of the weak 
mean magnetic field results in the small-scale magnetic 
buoyancy effect. We also investigate the anisotropic tur- 
bulent magnetic diffusion of the mean magnetic field in 
the turbulent convection. 

This paper is organized as follows. In section 2 we 
formulate the governing equations, the assumptions, the 
procedure of the derivations. In section 3 we consider 
the axisymmetric afl dynamo problem and determine the 
nonlinear drift velocities of the mean magnetic field and 
nonlinear turbulent magnetic diffusion in a turbulent con- 
vection. In section 4 we discuss the small-scale magnetic 
buoyancy and magnetic pumping effects and made esti- 
mates for the solar and stellar turbulent convection. Fi- 
nally, we draw conclusions in section 5. In Appendixes A 
we perform a detailed derivation of the nonlinear drift 
velocities of the mean magnetic field and nonlinear tur- 
bulent magnetic diffusion in a turbulent convection. 



2. The governing equations 

In this study we investigate the small-scale magnetic 
buoyancy and magnetic pumping effects in a turbulent 
convection. These phenomena are determined by the 
nonlinear drift velocities in the nonlinear electromotive 
force. In order to derive the nonlinear electromotive 
force in the turbulent convection we use a mean field 
approach in which the magnetic and velocity fields, and 
entropy are decomposed into the mean and fluctuating 
parts, where the fluctuating parts have zero mean val- 
ues. We assume that there exists a separation of scales, 
i.e., the maximum scale of turbulent motions Iq is much 
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smaller then the characteristic scale L of inhomogeneities 
of the mean fields. We adopt here a procedure of the 
derivation of the nonlinear electromotive force which was 
applied previously by Rogachevskii and Klccorin (2004) 
for the hydrodynamic incompressible turbulence. Let us 
outline here the procedure of the derivation of the non- 
linear electromotive force for the turbulent convection 
(for details, see also Appendix A). We consider a nonro- 
tating turbulent convection with large Rayleigh numbers 
and large hydrodynamic and magnetic Reynolds num- 
bers. The equations for fluctuations of the fluid velocity, 
entropy and the magnetic field are given by 
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where we use new variables (v, s, H) for fluctuating fields 
v = y/pou' and s = Jp~a~s', and also for the mean field 
H = B/(/iy/po). Here B is the mean magnetic field, po 
is the fluid density, /i is the magnetic permeability of the 
fluid, e is the vertical unit vector, fi 2 . = — g-VS 1 is the 
Brunt- Vaisala, frequency, S is the mean entropy, g is the 
acceleration of gravity, u', b and s' are fluctuations of 
velocity, magnetic field and entropy (for simplicity of no- 
tations we omitted prime in b because we did not use new 
variables for magnetic fluctuations), v , b N and s N are 
the nonlinear terms which include the molecular viscous 
and diffusion terms, p = p' + ^/po (H-b) are the fluctua- 
tions of total pressure, p' are the fluctuations of fluid pres- 
sure. Equations C3>-© f° r fluctuations of fluid velocity, 
entropy and magnetic field are written in the anelastic 
approximation, which is a combination of the Boussi- 
nesq approximation and the condition V- (po u ') = 0. 
The equation, V- u' = A p (u'-e), in the new variables 
reads: V-v = (A p /2)(v-e), where Vpo/po = — A^e. 
The quantities with the subscript "0" correspond to the 
hydrostatic nearly isentropic basic reference state, i.e., 
VP = pog and g-[(7^o) _1 VP - Po^Po] » 0, where 
7 is the specific heats ratio and Po is the fluid pressure in 
the basic reference state. The turbulent convection is re- 
garded as a small deviation from a well-mixed adiabatic 
reference state. 

Using equations CEJ-© written in a Fourier space we 
derive equations for the two-point second-order corre- 
lation functions of the velocity fluctuations (viVj), the 
magnetic fluctuations (bibj), the entropy fluctuations 
(ss), the cross- helicity (biVj), the turbulent heat flux 
(s Vi) and (s bi). The equations for these correlation func- 



tions are given by equations HA.5fl - ljA.10p m Appendix 
A. We split the tensor of magnetic fluctuations into non- 
helical, = (bibj), and helical, h\- , parts. The he- 

lical part h\j depends on the magnetic helicity and is 
determined by a dynamic equation which follows from 
the magnetic helicity conservation arguments (see be- 
low). We also split all second-order correlation functions, 
A/( /7 ), into symmetric and antisymmetric parts with re- 



hsj\ where 



spect to the wave vector k, e.g., 
the tensors h$ = [hij(k)+hij(— k)]/2 describes the sym- 



metric part of the tensor and h 



(a) 



Mk) 



-k)]/2 



determines the antisymmetric part of the tensor. 

The second-moment equations include the first-order 
spatial differential operators N applied to the third- 
order moments M( IIT \ A problem arises how to close 
the system, i.e., how to express the set of the third- 
order terms NM^ 111 ^ through the lower moments M^ 11 ^ 
(see, e.g., Orszag 1970, Monin and Yaglom 1975, Mc- 
Comb 1990). We use the spectral r approximation 
which postulates that the deviations of the third-moment 
terms, AfM^ 111 ^ (k), from the contributions to these 
terms afforded by the background turbulent convection, 
A/"M( /77 '°)(k), are expressed through the similar devia- 
tions of the second moments, M (II *>(k) - M (II '°^(k): 



AfM( IU \k) - AfM( nT > Xk) 



1 



[M< J/ )(k) 



r(k) 
-M( JI '°)(k)] 



(4) 



(see, e.g., Orszag 1970, Pouquet et al. 1976, Kleeorin et 
al. 1990, Kleeorin et al. 1996, Blackman and Field 2002, 
Kleeorin and Rogachevskii 2003, Rogachevskii and Klee- 
orin 2004, Brandenburg et al. 2004, Brandenburg and 
Subramanian 2005b, Kleeorin and Rogachevskii 2006), 
where r(k) is the scale-dependent relaxation time, which 
can be identified with the correlation time of the turbu- 
lent velocity field. In the background turbulent convec- 
tion, the mean magnetic field is zero. 

We apply the spectral r approximation only for the 
nonhelical part of the tensor of magnetic fluctuations. 

The helical part depends on the magnetic helicity, 
and it is determined by the dynamic equation which fol- 
lows from the magnetic helicity conservation arguments 
(see, e.g., Kleeorin and Ruzmaikin 1982, Gruzinov and 
Diamond 1994, Kleeorin et al. 1995, Gruzinov and Di- 
amond 1996, Kleeorin and Rogachevskii 1999, Kleeorin 
et al. 2000, Blackman and Field 2000, Kleeorin et al. 

2002, Blackman and Brandenburg 2002, Kleeorin et al. 

2003, Brandenburg and Subramanian 2005a, Zhang et al. 
2006, and references therein). The characteristic time of 
evolution of the nonhelical part of the tensor hij is of 
the order of the turbulent time tq — Io/uq, while the 
relaxation time of the helical part of the tensor h\- is 
of the order of To Rm, where Rm = IqUq/ti is the mag- 
netic Reynolds number (which is very large), uq is the 
characteristic turbulent velocity in the maximum scale 
of turbulent motions Iq and r\ is the magnetic diffusivity 
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due to electrical conductivity of the fluid. In this study 
we consider an intermediate nonlinearity which implies 
that the mean magnetic field is not enough strong in or- 
der to affect the correlation time of turbulent velocity 
field. The theory for a very strong mean magnetic field 
can be corrected after taking into account a dependence 
of the correlation time of the turbulent velocity field on 
the mean magnetic field. 

We assume also that the characteristic time of vari- 
ation of the mean magnetic field B is substantially 
larger than the correlation time r(fc) for all turbulence 
scales. This allows us to get a stationary solution 
for the equations for the second-order moments, M^ IT \ 
For the integration in k space of the second-order mo- 
ments we have to specify a model for the background 
turbulent convection which is determined by equations 
(|A.32jl - HA.34jl in Appendix A. This model takes into 
account the inhomogeneity of the turbulence described 
by the two parameters: A^ u) = V l (u /2 )(°V(u' 2 ) (0) and 



A 



Vi(b 2 )( )/{b 2 )(°). 



This model includes also 
(F) _ 



the inhomogeneity of the turbulent heat flux, A 



/<|u'| 



«'\(°) 



and the inhomogeneity of the 



fluid density described by the parameter A p . The 
quantities with the superscript (0) correspond to the 
background turbulent convection with B = 0. Us- 
ing the solution of the derived second-moment equa- 
tions, we determine the nonlinear electromotive force, 
£j = £imn J (b n v m )k dk, in the turbulent convection (see 
Appendix A), where Sijk is the fully antisymmetric Levi- 
Civita tensor. This allowed us to determine the nonlinear 
drift velocities of the mean magnetic field and nonlinear 
turbulent magnetic diffusion, and to study the small-scale 
magnetic buoyancy and magnetic pumping effects in the 
turbulent convection. 



3. The axisymmetric dynamo 

Let us consider the axisymmetric afl dynamo problem. 
The mean magnetic field in the local system of coordinate 
is B = B(x,z)e y + Vx [A(x, z)e y ], where B(x,z) and 
A(x, z) are determined by the dimensionless equations 
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dB 
~dt 



a(B)B + r/W(B 
-(V A (B)-V)A, 
D[V(ffl)xVi] 9 
-Vb(B)B] , 



d 2 A 
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d 2 A 
dx 2 



V-[f) B {B)VB 
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dn determine the differential rotation, D is the dynamo 
number (see below), a(B) is the total (hydrodynamic + 
magnetic) nonlinear alpha effect (see, e.g., Kleeorin et 
al. 2000, Rogachevskii and Kleeorin 2000, and references 
therein), r) is the diagonal tensor with the components 
^(2> x ) (B) f the nonlinear turbulent magnetic diffusion 
of toroidal field, rj^' x \B) are the nonlinear turbulent 



magnetic diffusion coefficients of the poloidal magnetic 
field, V^B) and V^(B) are the nonlinear drift velocities 
(see below). The axis z of the local system of coordinate 
is directed opposite to the gravity acceleration g and the 
axis x is in meridional plane and directed to the equator, 
so that the spherical coordinates (r, 8, <p) translate to the 
local system of coordinate (z,x, y). 

We adopt here the dimensionless form of the mean 
dynamo equations; in particular, length is measured in 
units of L, time is measured in units of the turbulent 
magnetic diffusion time L 2 /rj T and B is measured in 
units of the equipartition energy B cq = ■JfJipo Uq, a is 
measured in units of a* (the maximum value of the hy- 
drodynamic part of the a effect), the nonlinear turbu- 
lent magnetic diffusion coefficients are measured in units 
of rj T = IqUq/3, the nonlinear drift velocities ~Va,b(B) 
are measured in the units of r) T /L 1 the differential ro- 
tation Sfl is measured in units of <5f2* and the dimen- 
sionless parameters A^ u \ A^ b \ A p and A^ F ^ are mea- 
sured in the units of L~ x . We define R a = La*/ri T , 
= r (d(#n*)/dr) L 2 /rj T1 and the dynamo number 

D = RujRa- 

The derivation of equation for the nonlinear electromo- 
tive force allows us to determine the nonlinear turbulent 
magnetic diffusion coefficients and the nonlinear drift ve- 
locities of the mean magnetic field, which are given by 



r(«) 



a*^ 2) (B) 
a*V { A F B x) (B) 



V A , B (B) = V^ B (B) + a*V^(B) 



Here the superscript (v) corresponds to the contributions 
from the purely hydrodynamic turbulence and the su- 
perscript (F) corresponds to the contributions from the 
turbulent heat flux. These contributions are given by 
equations l|A.35() - l|A.47D in Appendix A. The parameter 
a» which is determined by the budget equation for the 
total energy, is given by 



1 



(8) 



where U is the mean velocity and v T is the turbulent 
viscosity. 

The asymptotic formulae for the nonlinear turbulent 
magnetic diffusion coefficients and the nonlinear drift ve- 
locities for the weak mean magnetic fields, B -C B oq /A, 
are given by 



vr'(B) = v, 



V ( B Z) (B) = 
1 + 0.1a* , 
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FIG. 1: Nonlinear turbulent magnetic diffusion coefficients 
(a) rj A and (b) r\ B in the vertical (solid) and horizontal 
(dashed) directions in a turbulent convection with a, = 0.8. 
The thin curves in figure lb correspond to e = and thick 
curves to e = 1. 



where we neglect the terms ~ 0(/3 2 ). Here (3 = 
V8B/B cq and the parameter e = (b 2 )(°7(v 2 ) (0) . When 
B 3> -Boq/4 the nonlinear turbulent magnetic diffusion 
coefficients and the nonlinear drift velocities are given by 
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FIG. 2: Vertical nonlinear drift velocities (a) V^ z ' and (b) 



Vg in a turbulent convection with a, = 0.8 (solid) and in a 
nonconvective turbulence, a* = (dashed) for A^ u ' = Ai F ' = 



A 



(S) 



A 



1. The thin curves correspond to e = 



and thick curves to e 



where we neglect the terms ~ 0(/3~ 2 ). 

The nonlinear turbulent magnetic diffusion coefficients 
^(z^z) Q £ ^ ue poioidal and toroidal components of the 
mean magnetic field in the vertical (along the z-axis) 
and horizontal (along the x-axis) directions are shown in 
figure 1 for the turbulent convection with a* = 0.8. The 
magnetic fluctuations due to the small-scale dynamo (de- 
scribed by the parameter e) increase the turbulent mag- 
netic diffusion of the toroidal mean magnetic field (see 
figure lb), and they do not affect the turbulent magnetic 
diffusion of the poloidal field. Note also that the nonlin- 
ear turbulent magnetic diffusion in a turbulent convec- 
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eg 
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FIG. 3: Vertical nonlinear drift velocities (V^ 
V^) in a turbulent convection with a* = 0.8 for A^ 



A 



(B) 



0, and for different depths h of the convective zone 



(from the solar surface) : (a) h = 1.7xl0 7 cm; (b) h = 3.7xl0 7 
cm; (c) h = 1.9 x 10 10 cm. The thin curves correspond to e = 
and thick curves to e = 1. 



tion is anisotropic even for a weak mean magnetic field. 
In particular, it is enhanced in the vertical (radial) direc- 
tion. 

The vertical nonlinear drift velocities of poloidal and 
toroidal components of the mean magnetic field in the 
turbulent convection (a* = 0.8) and in the nonconvec- 
tive turbulence (a* = 0) are shown in figure 2. The 
turbulent convection enhances the nonlinear drift veloc- 
ities of the mean magnetic field in comparison with the 
case of a purely hydrodynamic turbulence (sec figure 2). 
In the next section we discuss the nonlinear drift veloc- 
ities of the mean magnetic field in the solar convective 
zone which cause the small-scale magnetic buoyancy and 
magnetic pumping effects. 
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FIG. 4: Vertical nonlinear drift velocity V B ^ (figure 4a) of 
toroidal magnetic field in the overshoot layer with a„ = 0.8 



for A 



= Af > = 20, A p = A 



(B) 



1, A 



(b) 



A'"' - A fl 



The thick curve corresponds to e = 1, the thin solid curve 
corresponds to e = 0.9 and thin dashed curve corresponds to 
e = 0.5. Horizontal nonlinear effective drift velocities V^ B 
(figure 4b) of toroidal (solid) and poloidal (dashed) magnetic 
fields in the turbulent convection with a„ = 0.8 for = 
A^ = A^ F ' = A^ s ' = 1. The thin curves correspond to e = 
and thick curves to e = 1. 



4. Discussion 

Let us discuss the small-scale magnetic buoyancy and 
magnetic pumping effects. In figure 3 the vertical non- 
linear drift velocities of the toroidal and poloidal mag- 
netic fields are plotted for different depths h of the so- 
lar convective zone (measured from the solar surface): 



h = 1.7 x 10 7 cm (figure 3a); h = 3.7 x 10 7 cm (fig- 
ure 3b), and h = 1.9 x 10 10 cm (figure 3c). In order to 
estimate the governing parameters we use the models of 
the solar convective zone (see, e.g., Spruit 1974, Baker 
and Temesvary 1966). More modern treatments make 
little difference to these estimates. 

In particular, in the upper part of the solar convective 
zone, say at the depth h* ~ 1.7 x 10 7 cm, the parameters 
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are as follows: the characteristic turbulent velocity uq ~ 
2.2x 10 5 cm s _1 ; the maximum scale of turbulent motions 
l ~ 3.3 x 10 7 cm; the fluid density p ~ 4.6 x KT 7 g 
cm -3 ; the turbulent magnetic diffusion ?7 T ~ 2.4 x 10 12 
cm 2 s" 1 ; the density stratification scale A" 1 ~ 10 s cm 
and the characteristic scale of the inhomogeneity of the 
turbulent magnetic diffusion A" 1 = | V r rj T /-q T | _1 ~ 10 7 
cm. 

At the depth h* ~ 3.7 x 10 7 cm, the parameters are 
u ~ 1.5 x 10 5 cm s _1 ; l ~ 4.5 x 10 7 cm; p ~ 8.3 x 10~ 7 
g cm~ 3 ; 7y T ~ 2.3 x 10 12 cm 2 s _1 ; A" 1 ~ 4 x 10 7 cm 
and A" 1 ~ 2.2 x 10 8 cm. 

At the bottom of the solar convective zone, say at the 
depth h* ~ 1.9 x 10 10 cm, the parameters are u ~ 2x 10 3 
cm s _1 ; l ~ 8.1 x 10 9 cm; p ~ 2.1 x 10 _1 g cm" 3 ; 
rj T ~ 5.2 x 10 12 cm 2 s -1 ; A" 1 ~ 6.5 x 10 9 cm and 
A" 1 ~ 8 x 10 10 cm. 

Figure 3 demonstrates that only near the solar surface 
the radial drift velocity for a weak mean magnetic field 
is directed upward to the surface of the Sun. This causes 
to the small-scale magnetic buoyancy effect. However, in 
the main part of the solar convective zone the radial non- 
linear drift velocities of the toroidal and poloidal mean 
magnetic fields are directed downward. This results in 
the small-scale magnetic pumping effect. These phenom- 
ena are determined by the nonlinear drift velocities in the 
nonlinear electromotive force, and they are different from 
the large-scale magnetic buoyancy and magnetic pump- 
ing effects. The large-scale phenomena are caused by the 
effect of the mean magnetic field on the large-scale den- 
sity stratified fluid flow. These large-scale phenomena 
are stronger than the small-scale magnetic buoyancy and 
magnetic pumping effects when the mean magnetic field 
is larger than the equipartition field. In particular, the 
ratio of the velocities which correspond to the large-scale 
and small-scale effects, is of the order of (B/B cq ) 2 . 

In figure 4 the vertical (figure 4a) and horizontal (fig- 
ure 4b) nonlinear drift velocities of the toroidal mean 
magnetic field are plotted for the overshoot layer lo- 
cated at the bottom of the solar convective zone. In 
these layer the turbulence and the turbulent heat flux 
are strongly inhomogeneous. The drift velocities in fig- 
ures 2-4 are measured in the units of r\ T A p . Here we 

assumed that A* — A^ — A p , which implies that 

A| 6) = Vi(pa(u' 2 )W)/(p (u' 2 )W). Figure 4a demon- 
strates that the vertical nonlinear drift velocity of the 
toroidal mean magnetic field depends strongly on the 
level of the magnetic fluctuations caused by the small- 
scale dynamo (described by the parameter e). If there is 
a small deviation from e = 1 (the equipartition between 
the kinetic and magnetic turbulent energies) there is only 
the magnetic pumping effect in the overshoot layer. On 
the other hand, the horizontal nonlinear drift velocity of 
the toroidal mean magnetic field in the overshoot layer 
is negative, i.e., it is directed to the solar polar regions 
(see figure 4b). 

The magnetic pumping in three-dimensional compress- 



ible rotating magnetoconvection has been studied by Os- 
sendrijver et al. (2002) in direct numerical simulations 
(see also review by Ossendrijver 2003). The resulting 
pumping effects are isolated by calculating of the effective 
drift velocities in turbulent convection. The pumping ef- 
fects act differently on different components of the mean 
magnetic field (Ossendrijver et al. 2002). This result 
is in a good agreement with our results [see figure 2 and 
equations ©, @, fO?) . (Ogfr ]. The pumping ef- 
fect in the vertical direction is found to be equivalent to 
a predominating downward advection with a maximum 
drift velocity of the order of 10% of the turbulent velocity 
(Ossendrijver et al. 2002). This is in agreement with our 
theoretical findings (see, e.g., figures 2, 3 and 4a). Note 
that the effective drift velocity due to the inhomogeneity 
of the fluid density (see Kichatinov and Riidiger 1992) 
also causes a predominating downward drift of the mean 
magnetic field. 

The small-scale magnetic pumping and buoyancy ef- 
fects were investigated in the present study for large hy- 
drodynamic and magnetic Reynolds numbers using the 
spectral r approximation (the third-order closure pro- 
cedure). Previous analytical studies of the small-scale 
magnetic pumping and buoyancy effects (see Kichati- 
nov 1991, Kichatinov and Riidiger 1992, Kichatinov and 
Pipin 1993) were performed using the second order cor- 
relation approximation (SOCA). This approximation is 
valid for small hydrodynamic Reynolds numbers. In- 
deed, even in a highly conductivity limit (large magnetic 
Reynolds numbers) SOCA is valid only for small Strouhal 
numbers, while for large hydrodynamic Reynolds num- 
bers (fully developed turbulence) the Strouhal number 
is 1. In the present study we take into account the in- 
homogeneity of the fluid density assuming that (pu^u'A 
is weakly inhomogeneous (see equation ljA.320 ). This 
is in agreement with the models of the solar convec- 
tive zone (see, e.g., Spruit 1974, Baker and Temesvary 
1966). On the other hand, in studies by Kichatinov 
(1991) and Kichatinov and Riidiger (1992) it was as- 
sumed that (/9 2 u^-) is weakly inhomogeneous. Since the 
density in the solar convective zone varies over six or- 
ders of magnitude, the validity of the latter suggestion is 
questionable. 



5. Conclusions 

In summary, we study the nonlinear drift of the mean 
magnetic field in a turbulent convection. Three kinds of 
the inhomogeneities determine the nonlinear drift veloc- 
ities of the mean magnetic field: (i) the inhomogeneous 
turbulence; (ii) the nonuniform fluid density and (iii) the 
nonuniform turbulent heat flux. The inhomogeneous tur- 
bulence causes the well-known turbulent diamagnetic and 
paramagnetic velocities. The nonlinear drift velocities of 
the mean magnetic field result in the small-scale magnetic 
buoyancy and magnetic pumping effects in the turbulent 
convection. In the main part of the solar convective zone 
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the small-scale magnetic pumping effect dominates (i.e., 
the radial nonlinear drift velocity of the mean magnetic 
field is directed downward to the bottom of the convec- 
tive zone), while near the solar surface the small-scale 
magnetic buoyancy effect is important when the mean 
magnetic field is weak. These small-scale phenomena can 
be stronger than the large-scale magnetic pumping and 
magnetic buoyancy which are caused by the influence of 
the mean magnetic field on the stratified fluid flow. 

Acknowledgment 

An important part of this work was performed during 
our visit to the Isaac Newton Institute for Mathemati- 
cal Sciences (University of Cambridge) during the pro- 
gramme " Magnetohydrodynamics of Stellar Interiors" . 

References 

Baker N. and Temesvary, S., Tables of Convective Stellar 
Envelope Models, 1966 (NASA: New York). 

Blackman, E.G. and Brandenburg, A., Dynamic nonlin- 
earity in large scale dynamos with shear. Astrophys. J., 
2002, 579, 359-373. 

Blackman, E.G. and Field, G., Constraints on the mag- 
nitude of alpha in dynamo theory. Astrophys. J., 2000, 
534, 984-988. 

Blackman, E.G. and Field, G., New dynamical mean-field 
dynamo theory and closure approach. Phys. Rev. Lett., 
2002, 89, 265007 (1-4). 

Brandenburg, A., The case for a distributed solar dy- 
namo shaped by near-surface shear. Astrophys. J., 2005, 
625, 539-547. ' 

Brandenburg, A., Kapyla, P. and Mohammed, A., Non- 
Fickian diffusion and tau-approximation from numerical 
turbulence. Phys. Fluids, 2004, 16, 1020-1027. 

Brandenburg, A. and Subramanian, K., Astrophysical 
magnetic fields and nonlinear dynamo theory. Phys. 
Rept, 2005a, 417, 1-209. 

Brandenburg, A. and Subramanian, K., Minimal tau ap- 
proximation and simulations of the alpha effect. Astron. 
Astrophys., 2005b, 439, 835-843. 

Drobyshevski, E.M. and Yuferev, V.S., Topological 
pumping of magnetic flux by three-dimensional convec- 
tion. J. Fluid Mech., 1974, 65, 33-44. 

Galloway, D.J. and Proctor, M.R.E., The kinematics 
of hexagonal magnetoconvection. Geophys. Astrophys. 
Fluid Dyn., 1983, 24, 109-136. 

Gilman, P. A., Instability of magnctohydrostatic stellar 
interiors from magnetic buoyancy. Astrophys. J., 1970, 
162, 1019-1029. 



Gruzinov, A. and Diamond, P. H., Self-consistent theory 
of mean-field electrodynamics. Phys. Rev. Lett., 1994, 
72, 1651-1653. 

Gruzinov, A. and Diamond, P.H., Nonlinear mean-field 
electrodynamics of turbulent dynamos. Phys. Plasmas, 
1996, 3, 1853-1857. 

Kichatinov, L. L., Turbulent transport of magnetic fields 
in a highly conducting rotating fluid and the solar cycle. 
Astron. Astrophys., 1991, 243, 483-491. 

Kichatinov, L.L. and Pipin, V.V., Mean-field buoyancy. 
Astron. Astrophys., 1993, 274, 647-652. 

Kichatinov, L.L. and Rudiger, G., Magnetic-field advec- 
tion in inhomogeneous turbulence. Astron. Astrophys., 
1992, 260, 494-498. 

Kleeorin N., Kuzanyan, K., Moss D., Rogachevskii I., 
Sokoloff D. and Zhang, FL, Magnetic helicity evolution 
during the solar activity cycle: observations and dynamo 
theory. Astron. Astrophys., 2003, 409, 1097-1105. 

Kleeorin, N., Mond, M. and Rogachevskii, I., Magneto- 
hydrodynamic turbulence in the solar convective zone as 
a source of oscillations and sunspots formation. Astron. 
Astrophys., 1996, 307, 293-309. 

Kleeorin N., Moss D., Rogachevskii I. and Sokoloff D., 
Helicity balance and steady-state strength of dynamo 
generated galactic magnetic field. Astron. Astrophys., 
2000, 361, L5-L8. 

Kleeorin N., Moss D., Rogachevskii I. and Sokoloff D., 
The role of magnetic helicity transport in nonlinear galac- 
tic dynamos. Astron. Astrophys., 2002, 387, 453-462. 

Kleeorin, N. and Rogachevskii, I., Magnetic helicity ten- 
sor for an anisotropic turbulence. Phys. Rev. E, 1999, 
59, 6724-6729. 

Kleeorin, N. and Rogachevskii, I., Effect of rotation on 
a developed turbulent stratified convection: the hydro- 
dynamic helicity, the alpha-effect and the effective drift 
velocity. Phys. Rev. E, 2003, 67, 026321 (1-19). 

Kleeorin, N. and Rogachevskii, I., Effect of heat flux on 
differential rotation in turbulent convection. Phys. Rev. 
E, 2006, 73, 046303 (1-9). 

Kleeorin, N., Rogachevskii, I. and Ruzmaikin, A., Mag- 
netic force reversal and instability in a plasma with ad- 
vanced magnetohydrodynamic turbulence. Sov. Phys. 
JETP, 1990, 70, 878-883. Translation from Zh. Eksp. 
Teor. Fix., 1990, 97, 1555-1565. 

Kleeorin, N., Rogachevskii, I. and Ruzmaikin, A., Mag- 
nitude of dynamo - generated magnetic field in solar - 
type convective zones. Astron. Astrophys., 1995, 297, 
159-167. 



9 



Kleeorin, N. and Ruzmaikin, A., Dynamics of the av- 
eraged turbulent helicity in a magnetic field. Magne- 
tohydrodynamics, 1982, 18, 116-122. Translation from 
Magnitnaya Gidrodinamika, 1982, 2, 17-24. 

Krause, F. and Radlcr, K.-H., Mean-Field Magnetohy- 
drodynamics and Dynamo Theory, 1980 (Pergamon: Ox- 
ford). 

Kulsrud, R., A critical review of galactic dynamos. 
Annu. Rev. Astron. Astrophys., 1999, 37, 37-64. 

McComb, W.D., The Physics of Fluid Turbulence, 1990 
(Clarendon: Oxford). 

Moffatt, H.K., Magnetic Field Generation in Electrically 
Conducting Fluids, 1978 (Cambridge University Press: 
New York). 

Monin, A.S. and Yaglom, A.M., Statistical Fluid Me- 
chanics, 1975, Vol. 2 (MIT Press: Cambridge, Mas- 
sachusetts) . 

Moreno-Insertis, F., Schussler, M. and Ferriz-Mas, A., 
Enhanced inertia of thin magnetic flux tubes. Astron. 
Astrophys., 1996, 312, 317-326. 

Orszag, S.A., Analytical theories of turbulence. J. Fluid 
Mech., 1970, 41, 363-386. 

Ossendrijver, M., Stix, M., Brandenburg, A., Riidigcr, 
C, Magnetoconvection and dynamo coefficients. II. 
Field-direction dependent pumping of magnetic field. 
Astron. Astrophys., 2002, 394, 735-745. 

Ossendrijver, M., The solar dynamo. Astron. Astrophys. 
Rev., 2003, 11, 287-367. 

Parker, E., The formation of sunspots from the solar 
toroidal field. Astrophys. J., 1955, 121, 491-507. 

Parker, E., The dynamical state of the interstellar gas 
and fluids. Astrophys. J., 1966, 145, 811-833. 

Parker, E., Gosmical Magnetic Fields, 1979 (Oxford Uni- 
versity Press: New York). 

Pouquet, A., Frisch, U. and Leorat, J., Strong MHD tur- 
bulence and the nonlinear dynamo effect. J. Fluid Mech., 
1976, 77, 321-354. 

Priest, E.R., Solar Magnetohydrodynamics, 1982 (D. Rei- 
del Publ. Co.: Dordrecht). 

Radler, K.-H., Kleeorin, N. and Rogachcvskii, I., The 
mean electromotive force for MHD turbulence: the case 
of a weak mean magnetic field and slow rotation. Geo- 
phys. Astrophys. Fluid Dynam., 2003, 97, 249-274. 

Roberts, P.H. and Soward, A.M., A unified approach to 
mean field electrodynamics. Astron. Nachr., 1975, 296, 
49-64. 



Roberts, P.H. and Soward, A.M., Dynamo theory. Annu. 
Rev. Fluid Mech., 1992, 24, 459-512. 

Rogachevskii, I. and Kleeorin, N., Electromotive force 
for an anisotropic turbulence: intermediate nonlinearity 
Phys. Rev. E, 2000, 61, 5202-5210. 

Rogachevskii, I. and Kleeorin, N., Nonlinear theory of a 
"shear-current" effect and mean-field magnetic dynamos. 
Phys. Rev. E, 2004, 70, 046310 (1-15). 

Ruzmaikin, A. A., Shukurov, A.M. and Sokoloff, D.D., 
Magnetic Fields of Galaxies, 1988 (Kluwer Academic: 
Dordrecht). 

Schussler, M., Caligari, P., Ferriz-Mas, A. and Moreno- 
Insertis, F., Instability and eruption of magnetic flux 
tubes in the solar convection zone. Astron. Astrophys., 
1994, 281, L69-L72. 

Spiegel, E.A. and Weiss, N.O., Magnetic activity and 
variations in solar luminosity. Nature, 1980, 287, 616- 
617. 

Spruit, H.C., A model of solar convective zone. Solar 
Phys., 1974, 34, 277-290. 

Spruit, H. O, Motion of magnetic flux tubes in the solar 
convection zone and chromosphere. Astron. Astrophys., 
1981, 98, 155-160. 

Spruit, H.C. and van Ballegooijen, A. A., Stability of 
toroidal flux tubes in stars. Astron. Astrophys., 1982, 
106, 58-66. 

Stix, M., The Sun: An Introduction, 1989 (Springer: 
Berlin and Heidelberg). 

Thomas, J.H., Weiss, N.O., Tobias, S.M. and Brummell, 
N.H., Downward pumping of magnetic flux as the cause 
of filamentary structures in sunspot penumbrae. Nature, 
2002, 420, 390-393. 

Tobias, S.M., Brummell, N.H., Clunc T.L. and Toomrc, 
J., Transport and storage of magnetic field by overshoot- 
ing turbulent compressible convection. Astrophys. J., 
2001, 549, 1183-1203. 

Tobias, S.M. and Hughes, D.W., The influence of veloc- 
ity shear on magnetic buoyancy instability in the solar 
tachocline. Astrophys. J., 2004, 603, 785-802. 

Vainshtein, S.I. and Kichatinov, L.L., The macroscopic 
magnetohydrodynamics of inhomogencously turbulent 
cosmic plasmas. Geophys. Astrophys. Fluid Dynam., 
1983, 24, 273-298. 

Weiss, N.O., The expulsion of magnetic flux by eddies. 
Proc. R. Soc. Lond., 1966, A293, 310-328. 

Zeldovich, Ya.B., The magnetic field in the two- 
dimensional motion of a conducting turbulent fluid. Sou. 



10 



Phys. JETP, 1957, 4, 460-462. Translation from Zh. 
Eksp. Teor. Fiz., 1956, 31, 154-156. 

Zeldovich, Ya.B., Ruzmaikin, A. A. and Sokoloff, D. 
D., Magnetic Fields in Astrophysics, 1983 (Gordon and 
Breach: New York). 

Zhang, H., Sokoloff, D., Rogachevskii, I., Moss, D., Lam- 
burt, V., Kuzanyan, K. and Kleeorin N, The Radial Dis- 
tribution of Magnetic Helicity in Solar Convective Zone: 
Observations and Dynamo Theory. Mon. Not. R. As- 
tron. Soc, 2006, 365, 276-286. 



APPENDIX A: THE NONLINEAR 
ELECTROMOTIVE FORCE IN TURBULENT 
CONVECTION 

Let us derive equations for the second-order moments 
in a turbulent convection. For this purpose we rewrite 
equations Q)-© in a Fourier space. In particular, 



d^(k, t) 
At 



E lm S$(b; H) + g e„ P m (k) s(k, t) 



iA, 



+A P D imn (k) S%l(b; H) + —%gk m Pim{e) s(k, t) 



2fc 2 



(A.l) 



dbi(k,t) A, 



- -f Rimn S%l(v; H)+ik m S%}(v; H) 



-S c >(v;H)+b? , 



(A.2) 



where we multiply equation Q written in k-space by 
Pij (k) = Sij — kij in order to exclude the pressure term 
from the equation of motion, and 



slfi^A) = J a,(k-Q)A i (Q)dQ, 
S^(a; A) = (2P m (fc) - S m )Si c \a; A) + ik n S$(a; A) , 
^ (c) (a;A) = i /a p (k-Q)g p A 4 (Q)dQ , 



Eij = S i: j - {\K p /k 2 ){k i e j - %(k-e)) , 

P*imn — &pPip(ty3mn ~t~ ^-p^mpSin (l/2)Cn^im j 
Pimn — &m,0~in &n$im j 

Pij (k) — Sij — k^ , Pij (e) = 5ij — eij , Sij is the Kronecker 
tensor, kij = kikj/k 2 and = CiCj. Here we neglect 
terms ~ 0(A 2 ). We use the mean- field approach, and the 
two-point correlation function of the velocity fluctuations 
is given by 

(vi(x)Vj(y)) = J / y (k, R) exp (ik-r) dk , 

where hereafter we omit argument t in the correlation 
functions, fij(k, R) = L(vi;vj), and 



L(a; c) = / (a(t, k + K/2)c(t, -k + K/2)) 

xexp(iK-R)dK , (A.3) 



(see Roberts and Soward 1975). Here R — (x + 
y)/2, r = x — y. Note that R and K correspond to 
the large scales, and r and k to the small ones. 

Using equations (|A.1(1 - I|A.2|) we derive equations 
for the following correlation functions: 



f i3 (k) = L{vi\Vj) , hij(k) = L{bi\bj) 
Fi(k) = L(s;vi) , ffijOO = L(bi,Vj) , 
Gi(k) = L(s;bi) , 0(k) =L(s;s) . 

These equations are given by 



dt 
dhjjjk) 

dt 
% 3 ( k ) 

dt 

dF t (k) 
dt 

dG t (k) 
dt 

de (k) 

dt 



(A.4) 

(A.5) 
(A.6) 



= :;k.TI !•!.,, ■ • XI , , 

= iik- ii • /;; • A'// (/ , 

= i(k-H)[/ y -(k)-Mk)] 

+ge n P jn {k)Gi(-k) + If. +N 9l] ,(A.7) 

= -i(k-H)G i (k)+ ff e n P in (fc)e(k) 

+/f+A r F i , (A.8) 

= -i(k-H)Fi(k) + if +NGi , (A.9) 
O 2 

= -F z (k) +A"6 , (A. 10) 



where hereafter we also omit argument R in the correla- 
tion functions. Here $y (k) = (k) — gji{— k), and 



i£=i£.(k) + i*(-k) 



JV^-(k) + Mii^(k) , j*(k) = N? n9jn (-k), 
</ BJ -(k) + J\#,Mk) + A^Gi(-k) , 
iV/ n G n (k) - M 4 0(k) , if = Af„F„(k) , 
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Mi 



d 



-Hi. 



2A; 2 

Apf i n (e)fc n ] , 



V - 9/9R and fly = Vjfl;, = ffe tl [P iri (fc)J , J -(k) + 
P,- n (*)Fi(-k)] +^4-, and A/%-, Nhij, Af gij , NF U 
NGi and A/ Q are the third-order moment terms appear- 
ing due to the nonlinear terms. The terms ~ Fi in the 
tensor J\fij can be considered as a stirring force for the 
turbulent convection. Note that a stirring force in the 
Navier-Stokes turbulence is an external parameter. 
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For the derivation of equations (|A.5(1 - I|A.10(1 we use an 
identity for the function (k, R) : 

Z y (k, R) = i J(k p + K p /2)H p (Q) exp(iK-R) 
x(«j(k + K/2-Q)^(-k + K/2))dKdQ . 
The identity reads 



Z«(k,R) ~ [i(k.H) + (l/2)(H-V)]/. y (k,R) 



l_ k S/ y (k) 



dk. 



(A.ll) 



(see Rogachevskii and Kleeorin 2004), and similarly for 
other second-order moments. We take into account that 
in equation (|A.7J) the terms with symmetric tensors with 
respect to the indexes "i" and "j" do not contribute to 
the nonlinear electromotive force. In equations i|A.5|) - 
we neglect the second-order and high-order spatial 
derivatives with respect to the large-scale variable R. 

Let us solve equations (|A.5jl - (|A.10|l neglecting the 
sources ... with the large-scale spatial deriva- 

tives. Then we take into account the terms with the 
large-scale spatial derivatives by perturbations. Thus, 
subtracting equations I|A.5|I - (|A.10J| written for back- 
ground turbulent convection (i.e., for B = 0) from those 
for B / 0, using the spectral r approximation [which 
is determined by equation @], neglecting the terms 
with the large-scale spatial derivatives, assuming that 
rjk 2 <C t _1 and vk 2 <C t -1 for the inertial range of turbu- 
lent fluid flow, and assuming that the characteristic time 
of variation of the mean magnetic field B is substantially 
larger than the correlation time r(k) for all turbulence 
scales, we arrive to the following steady-state solution of 
the obtained equations: 



4 0) (k)+ir(k.H)^(k) 
ftg>(k)-ir(k.H)*«(k) 
ir(k-H)[4-(k)-^(k)] 



hij(k) 

m) 

G t (k) 

e(k) 

where fy, hij, 

l( p I*,..., If and <Mk)=^(k) 



Tge n P jn (k)G t (k) , 
; (0) (k)-ir(k-H)G 2 (k) 
T 9 e n P in (k)[Q(k) - Q^(k)] 
ir(k-H) J F 4 (k) , 



e< '(k) 



(A.12) 
(A.13) 

(A.14) 

(A.15) 
(A.16) 
(A.17) 



9 are the solutions without the sources 
/., i£, ...,/f and $y (k) = gij(k) - foj(- k). Th e quan- 
tities with the superscript (0) in equations ljA.12(l - l|A.17|) 
correspond to the background turbulent convection. Here 
we take into account that for the background turbulent 
convection glf(k) = and Gf\k) = 0. 

Now we split all second-order correlation functions into 
symmetric and antisymmetric parts with respect to the 



wave vector k, i.e., 



flf +m\ where / 



t (a) 



[/y(k) + &(-k)]/2 and fjj> = [/^(k) - M-k)]/2. 
Thus, equations l|A.12|) - IIA.17|l yield 



/J } (k) 



^(k) 



1 



-2^rge n P in (k)F^(k)] 
H(l + 



g s) ( k ) 

(A.18) 



1 



1 



;^4° s) (k) 



£«(k) 



2ip [ 

+^Tge n P m {k)F^{k)] 
ir(k-H) (0fl) ( 0s) 
1 + 2^ [hj { ' ij 
+Tge n P m (k)F ( j s \k)] , 

,(0«) 



^ S \k) 



(k) 



l + V/2 ' 
Gr J (k) « -ir(k-H) J F;W(k), 

where ip(k) = 2(rk-H) 2 and we neglect terms 
in equations 1A.17I) . The correlation functions /, 



(A.19) 

(A.20) 
(A.21) 

(A.22) 
0(0?) 

(a) 



h) A ,g\j F> ' and G\ ' vanish if we neglect the large- 



Co) p(a) 
ij ' &ij i 

scale spatial derivatives, i.e., they are proportional to the 
first-order spatial derivatives. 

Now we take into account the large-scale spatial deriva- 
tives in equations i|A.5|) - (|A.10l) by perturbations. Their 
effects determine the following steady-state equations for 
the second moments: 



£ 7 o) (k) 



^(k) 



£§ ] (k) 



/J?° 3 (k) 

+Tl h - 



ir(k-H)*W(k) 



ir(k-H)$g>(k) 



= r[i(k-H)(/W(k)-/iW(k)) 
+ g e n P jn (k)G\ s \-k)+If j ] 
G\ s '(k) = -r[i(k.H)^(k)-7f], 



;(»)/ 



^ (o) (k) 



= ^(k)-r[i(k-H)Gr ; (k) 
-If] , 



(*), 



(A.23) 
(A.24) 

(A.25) 
(A.26) 

(A.27) 



where the second moments , ftjj gij , ... determine 



the effect of the large-scale derivatives and ^>[f (k) 



(k) — g^f(-k). The correlation functions of the 



background turbulent convection / <7 - (k), /i Vjj 10 '' (k) and 



,(0o) 



j^(° a ) ^ are determined by the inhomogeneities of turbu- 
lence, the fluid density and the turbulent heat flux [see 
equations (|A.32(l - i|A.34(l below]. Equations l|A.26|) and 
(lA~2TI) yield 



ff>(k) 



G, (s) (k) 



l + ip/2 



(0o) 



1 + VV2 



i(k-H)r/f 



[i(k.H)(^ (0a) +rJf) 



(A.28) 



12 



-If] 



(A. 29) where 



Our goal is to calculate the mean electromotive force 
£i(r = 0) = (1/2^/ pp) £jnm f im l(k) dk. Solution of sys- 
tem of equations (|A.23(l - i|A.25(l allow us to get the ex- 
pression for (k) which yields the mean electromotive 
force: 



& = 



TSi 



i(k.H)[/W-^)+r(J^ 



mn rtin/ 



+i 9 

1 mn 



rffe p P mf ,(fc) ^ 0Q) 
1 + ^/2 L V A ™ 



I G 1 

± n J 



dk 



where we use equations (IA.28|) and (|A.29|) 
tion iTOOll can be rewritten in the form: 



(A.30) 
Equa- 



& = 



1 + 2^ 



i(k.H){/i°„ a )-^) + 2r[(7V/ 



+N*)gpn ~ M n (F m - F^)}} + NtJpn 



Nl h pn + iTM n (k-B.)F„ 



i(k.H)F(°°)+(|iV/ q -N*)F q 



l + V/2 



dk . 



(A.31) 



For the integration in k-space in equation (|A.31(1 we spec- 
ify a model for the background turbulent convection (i.e., 
the turbulence with zero mean magnetic field, B = 0), 
which is determined by 



#(k) 



/4?(k) 



*f>Ck) 



/* W(k) 
h*W(k) 



^■(k) + ^(fc*Af 



(6), 



(A.32) 



(A.33) 



3F*W{k) ej 



(A.34) 



0(°)(k) = 2e*W(k), g\y(k) = and G\°>(k) = 0, 
where -Py(k) = Sj,j — fcy, fey = kikj/k 2 , W(k) = 
E(k)/8Trk 2 , r(fc) = 2r f(fc), B(jfc) = -df(k)/dk, 
f(k) — (k/k ) 1 ^ q , 1 < q < 3 is the exponent of 
the kinetic energy spectrum (e.g., q — 5/3 for Kol- 
mogorov spectrum), fco = l/'o and to = Iq/uq. Here 
A<* > = A« - 2A p e, and Af > = Af > 
imply that 



2Ape 4 . These 



A 



W _V^ (u' 2 )W) 



,m(o) 



pi(|u'| S ') (0) 



(») 



A 



A„ (i?) = 



V,(u /2 )(°) 

(U' 2 )(°) ' 

Vi(|uV)(°> 
(|u'|s') (0) 



A 



(6) 



Vi(b 2 )W 
(b 2 )(°) ' 



and r*f J (k)dk = FmBi , J7^(k)dk = (/,/3)%, 

/ hf) (k) dk = (hm/3) S i:j and / 0(°) (k) dk = 9*. 

After the integration in k space in equation (|A.31Jl we 
obtain the nonlinear electromotive force. This yields the 
nonlinear turbulent magnetic diffusion coefficients and 
the nonlinear drift velocities of the mean magnetic field in 
the axisymmetric case, which are given by equation Q, 
where the contributions from the purely hydrodynamic 
turbulence are given by 



V£>(B) 



= 4 X) (4B) + A£' (45) 
= 4 1) (4S) + 3(l-e) 



-— A 2 {1QB 2 

Z7T 



(A.35) 



4 1} (4B) 



(A.36) 



A( B ) r m 
.__^(2-3e)AW(4B) 

3(1 -e) 



v£>(5) 



2vr 



A 2 (16S 2 ) 



»jW(fl)(A 



(u) 



2 



2 p 



eA< 6 >) 



(A.37) 



(A.38) 



e4 1) (4S)-(5-6e)A 2 1J (4S) 

(A.39) 



+ ^^A 2 (16B 2 ) 



2?r 



and the contributions caused by the turbulent heat flux 
are 



2*i{A 2 - 3Ai + 3d} + 4* 2 {Ai - Ci} 



e- 2) (s) = 



-3* 3 {Ai + Ci} , (A.40) 

2*!{A! +Ci} + 4* 2 {Ci} 
2,^{A 1 -2C 1 ) , (A.41) 

(-6*1 + 4* 2 + 3* 3 ){^i + A 2 - d 
C 3 } + 2* 1 {C 3 } + 6* 3 {Ci} , 

(A.42) 
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(A.43) 
(A.44) 



4 z 



y 4 {A 1 +A 2 -C 1 -C 3 } 
-A('{4B) - Af\AB) +C[ 2 \AB) 



+cf (4B) 



(A.45) 



VP{B) = 



^ 4 {^i + A 2 - 5d - 5C 3 } 



+5c { 3 2) m 



(2), 



(A.46) 



*i{17Ai + 17^2 + VTCi - 7C 3 } 
* 2 {6Ai + 6A 2 + lOCi - 2C* 3 } - ^{^i 



-9A 2 + 15Ci - 3C 3 } - 44 2) (4£) 



44 2) (4S) 



(A.47) 



where = (VB 2 )/B 2 , the parameter e = 

(b 2 )(°) /(v 2 )^°\ the parameter a* is given by equation (JHJl, 



and 



*i{A} 
* 2 {A} 



* 3 {A} = -- 
18 



\[AX^ {AB) - X [2 \2B)] , 
16A (2) (4S) - 16A (2) (2S) 

A(4B 2 ) 

4A (2) (4B) -49A (2) (2B) 



* 4 {X} 



— X{AB Z 



3XW(4B) - —X(16B 2 ) . (A. 

27T 



Note that * 4 {Ai} = A{ 1} + (l/2)A\ L > . The functions 

4m HP), C&\/3) for n = 1; 2 and the functions ^ m (/5 2 ), 
C m ((3 2 ) are given in Rogachevskii and Kleeorin (2004, 
in Appendixes B, C and D). Asymptotic formulae for 
the nonlinear turbulent magnetic diffusion coefficients 
and the nonlinear drift velocities of the mean magnetic 
field in the axisymmetric case are given by equations @ 
and (UJ. 
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